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Résumé: Une comparaison est faite entre le modèle 
fluide et la simulation de monte Carlo d'une décharge 
radiofréquence dans l'argon L'accord dans les 
propriétés prédites par les deux modèles est examiné 
pour une pression de 1 torr et une tension de 300 v . Un 
accord quantitatif a été trouvé dans le comportement 
général de la décharge prédit par le modèle fluide et la 
simulation Monte Carlo à cette pression et cette tension. 


L. INTRODUCTION 


Les application industrielles des décharges 
électriques ont trouvé une large diffusion dans différents 
secteurs industriels . Elles ont envahi des domaines 


d'importance vitale, tels que l'industrie micro- 
électronique, en faisant fin aux techniques 
traditionnelles. 


L'utilisation des décharges luminescentes radio- 
fréquences (> 100 Khz) à basse ou moyenne pression est 
un moyen très efficace pour créer des plasmas à basse 
température (plasmas froids) . De tels plasmas réactifs 
sont utilisés dans la gravure et le dépôt des couches 
minces qui représentent deux étapes fondamentales dans 
la fabrication des circuits intégrés . 

La connaissance des mécanismes physiques et 
chimiques se déroulant dans le plasma des décharges 
radiofréquence est indispensable pour l'optimisation des 
réacteurs à plasmas . Une meilleure compréhension de 
ces mécanismes complexes ne peut être obtenue que par 
la modélisation . 

En effet, plusieurs modèles ont été développés dans 
ce contexte. Ils peuvent être classés en deux catégories : 
les modèles microscopiques ou cinétiques[1] et les 
modèles macroscopiques ou fluides [2] - [10]. 

Les modèles microscopiques (cinétiques) sont basés 
sur la résolution directe de l’équation de Boltzmann par 
des méthodes des différences finies ou des éléments finis 
ou indirectement par des simulations de Monte Carlo. 
Comme inconvénient, ces modèles sont beaucoup 
consommateurs de temps surtout à pressions élevées. 

Les modèles macroscopiques ou fluides, utilisés à des 
pressions élevées, font intervenir à la place de l’équation 
de Boltzmann les moments d’ordre zéro, un et deux de 
cette équation couplés à l’équation de Poison. Ces 
modèles sont rapides . 

Les modèles cinétiques donnent la fonction de 
distribution des vitesses ou des énergies des particules 
chargées . Par contre, les modèles macroscopiques 
donnent des valeurs moyennes des densités et des 
vitesses moyennes des particules chargées, en supposant 
une distribution maxwelliane des énergies des électrons . 
Dans ce travail, on présente une comparaison des 
prédictions du modèle fluide et celles d'une simulation 


de Monte Carlo d'une décharge luminescente 
radiofréquence dans l'Argon (gaz électropositif) à une 
pression de 1 torr. La formulation des modèles est 
décrite dans la section 2 . La section 3 contient les 
résultats et les prédictions de chaque modèle. 
Finalement, une conclusion est donnée dans la section 4. 


II. DESCRIPTION DES MODELES 


A. Simulation de Monte Carlo 


La simulation de Monte Carlo utilisée dans ce travail est 
décrite en détail dans la référence [1]. Brièvement, cette 
méthode consiste à simuler les collisions entre particules 
(électrons-atomes , ions-atomes) par un tirage de 
nombres aléatoires et retracer l'historique de chaque 
particule depuis sa création jusqu'à sa disparition. 
Les électrons commencent leur mouvement à partir de la 
cathode avec une énergie, en moyenne, égale à 4 eV. la 
direction de chaque électron, définie par l'angle axial 6 et 
l'angle radial ®, est donnée par: 
cosO=R, , p=27r-R; 
où R, et R, sont des nombres tirés de l'intervalle (0,1). 

Les vitesses initiales sont calculées à partir de 
l'énergie E et les angles 6 et o : 

2E 


Vo = SinOsinp ru (1) 
Vo = cos8./2E 


Le mouvement des électrons, sous l'action du champ 
axial E, , est déterminé par les lois de la mécanique 
classique : 

X=Xp+ Vo - At 


Y = Yo + Vyo At (2) 
Z= 20 + Vo At (At) 

et Vk = Vo 
MEN G) 
Ve = Vo At 


OÙ Xo, Yo, Zo €t X, y, Z sont les coordonnées des électrons 
aux instants t et t + At respectivement, Vs, Vo, Vo et 
V,, V,, V, sont les vitesses correspondantes et m. et e 
sont respectivement la masse et la charge de l'électron. 
Le champ axial E, est obtenu par la résolution de 
l'équation de Poisson : 
SP (4) 
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p étant la charge d'espace et & la permittivité du vide. 

Après l'expiration du temps At (pas de temps), la 
probabilité de collision est calculée par un tirage d'un 
nombre aléatoire R;. Si une collision a lieu, le type de 
collision (élastique, ionisation, excitation) est déterminée 
par un tirage d'un nombre aléatoire R4 .Dépendant du 
type de collision, les nouvelles énergie et directions sont 
calculées . Le mouvement d'électron, après le choc et 
pendant un pas de temps At, est décrit de nouveau par les 
lois de Newton. Cette procédure est répétée jusqu'à ce 
que cet électron entre en collision avec les parois. 

Dans le cas de collisions ionisantes, le transport des 
nouveaux électrons et les ions positifs correspondant est 
décrit, de la même façon que précédemment, par la 
méthode de Monte Carlo. Seules les collisions de type 
élastique et transfert de charge sont considérées dans la 
simulation du mouvement des ions positifs. 


B. Modèle fluide 


Le modèle fluide utilisé dans cette étude est 
unidimensionnel et similaire à celui utilisé par Gogolides 
dans le référence [8]. Les équations issues des trois 
premiers moments de l’équation de Boltzmann sont: 


— les équations de continuité des électrons et des 
ions positifs : 


ne À 
et + cu = k;Nne (5) 
on ôj 
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où ne , lp ; Je Et JP représentent les densités et les flux des 
électrons et des ions positifs. N étant la densité du gaz 
que l'on suppose rester constante. k; représente le 
coefficient d’ionisation. 

— les équations de transfert de quantité de mouvement 
des électrons et des ions positifs : 


; °n 
Je = De: _ — HeneË (7) 
; On ef 
ÿ =D, En +UpnpE (8) 
D. , D, , He et u, représentent respectivement le 


coefficient de diffusion et la mobilité des électrons et des 


ions positifs. ET est le champ que voient réellement les 
ions positifs , car le mouvement de ces derniers 
(particules lourdes) est en retard de phase par rapport au 


champ électrique appliqué E. Ce champ effectif E* est 
exprimé par l'équation [2] : 


ef 
= L vT(E _. É) (0) 
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où v =—— représente la fréquence effectif de 


Mimi 
transfert de quantité de mouvement. 
— l'équation d'énergie des électrons: 
Ô(n,E) ” êq 
t x 


= €Eje — ki Nn, (10) 


où € représente l'énergie moyenne des électrons. k, étant 


le coefficient de perte d’énergie des électrons. q 
représente le flux d'énergie électronique et peut être 
donné par la formule approchée suivante[8]: 
s dE 
ape 
Ces équations doivent être couplées à l'équation de 
Poisson: 


dv e 
4 (11) 
E = = (12) 


Le coefficient d’ionisation k; peut être exprimé en 
fonction de l’énergie moyenne des électrons par [8]: 
1 
k,=8,7x10 °(&—5,3)exp(-49/(&-5,3)2) (13) 
pour & > 5,3ev. k; est supposé nul pour € < 5,3ev. 


Le coefficient de perte d'énergie kL peut être aussi 
exprimé en fonction de l'énergie moyenne des électrons 
par [8]: 
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Nous considérons une décharge radiofréquence à 
électrodes parallèles qui s'étend de x=0 à x=L. On 
applique à l'électrode de gauche (x=0). un potentiel rf:: 
v(t)= v,, sin(® çt). L'électrode de droite est mise à 
la terre (x=L). Les densités des électrons, des ions 
positifs sont supposées nulles aux électrodes. L'énergie 
électronique est supposée égale Lev aux électrodes 


IIT. RESULTATS 


Dans cette section, nous présenterons les résultats 
des deux modèles (fluide et cinétique) d'une décharge rf 
dans l'argon (gaz électropositif) à une pression de 1 
Torr. Une tension de 300v de crête et 10 Mhz de 
fréquence est appliquée à l'électrode de gauche (x = 0). 
L'électrode de droite (x = L) est mise à la terre. L'espace 
entre électrodes est pris égal à 2 cm . L'émission 
électronique secondaire aux électrodes est supposée 
nulle (y=0). 

Les paramètres de transport de l'argon utilisé dans le 
modèle fluide (tirés de l'article de Gogolides[8]), sont : 


D, = 120x10* cm? .Torr/s ; He = 30x10% cm° .Torr / v.s 
D, = 40.0 em”. Torr /s et p, = 14.44x10° em°.Torr / v.s. 

Les sections efficaces d'ionisation , d'excitation et des 
collisions élastiques électron-atome utilisées dans la 
simulation de Monte Carlo sont calculées par le code 
BOLZIG [11]. 


La section efficace de transfert de charge ion-atome 
est prise constante et égale à 10° cm° [12]. Quant à la 
section efficace de collision élastique ion-atome est 
29.10 '° cm’ (calculée). 


Les figures 1(a) et 1(b) représentent respectivement 
les variations spatio-temporelles de la densité des 
électrons prédites par le modèle fluide et la simulation de 
Monte Carlo. Les deux modèles prédisent la même 
forme et la même allure des variations de la densité 
électronique. Le balayage des électrons vers le plasma se 
produit au mêmes moments du cycle rf dans les deux 
simulations. La densité électronique, dans le plasma, est 
de 6x10° cm” d'après le modèle fluide et de 4x10° cm” 
dans la simulation de Monte Carlo. Ce qui prouve qu'il y 
a un accord quantitatif et qualitatif entre les deux 
modèles dans la description des électrons. 
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FIG.1: Variations spatio-temporelles de la densité des électrons dans 
une décharge luminescente radiofréquence dans l'argon pour les 
conditions suivantes: p = 1 torr, V = 300 v, f= 10 Mhz, et d =2 

cm., a) simulation de M.C, b). modèle fluide 


Les figures 2(a) et 2(b) présentent les variations 
spatio-temporelles du champ électrique. Les deux 
modèles prédisent la même distribution du champ 
électrique que ce soit dans les gaines ou dans le plasma. 
Ils prédisent, aussi, le même signe et presque la même 
amplitude du champ électrique pendant le balayage des 
électrons du ou vers le plasma. Ce qui permet d'affirmer 
d'avantage l'accord , entre le modèle fluide et la 
simulation de Monte Carlo à ces conditions. 
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FIG.02: Variations spatio-temporelles du champ électrique dans 
une décharge luminescente radiofréquence dans l'argon pour les 
conditions suivantes: p = 1 torr, V = 300 v , f = 10 Mhz, et d = 
2 cm., a) simulation de M.C, b) modèle fluide 


Les deux modèles prédisent le même courant de 
déplacement (voir figures 3 (a) et 3 (b) ) puisqu'ils 
donnent la même distribution du champ électrique. En 
effet, presque la même valeur maximale du courant de 
déplacement dans les gaines (3,6 mA/cm° pour le 
modèle fluide et 3 mA.cm” pour la simulation de M.C.) 
est prédite par les deux modèles. Le courant de 
déplacement dans le plasma est quasiment nul dans les 
deux modèles (champ faible et quasiment constant dans 
le plasma). 
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FIG.04: Variations spatio-temporelles du courant électronique 
dans une décharge luminescente radiofréquence dans l'argon 
pour les conditions suivantes: p = 1 torr, V = 300 v, f= 10 Mhz, 
et d=2 cm. a) simulation de M.C, b) ) modèle fluide. 


Le profile du courant électronique présente la même 
allure dans les deux modèles (voir fig. 4 (a) et 4 (b) ). Il 
est important dans le plasma (densité électronique 
maximale dans le plasma). et faible ou même nul dans 
les gaines (électrons balayés vers le plasma). 
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FIG.04: Variations spatio-temporelles du courant électronique dans 
une décharge luminescente radiofréquence dans l'argon pour les 
conditions suivantes: p = 1 torr, V = 300 v , f = 10 Mhz, et d=2 

cm., a) simulation de M.C ,b) modèle fluide 


Quant au courant ionique, les deux modèles prédisent 
des valeurs relativement faibles dans les gaines(0.06 
mA/cm° dans la modélisation fluide et 0.04 mA/cm° dans 
la simulation de M.C.) et des valeurs plus faibles 
(presque nulles) dans le plasma. 
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FIG.05: Variations spatio-temporelles du courant ionique dans 

une décharge luminescente radiofréquence dans l'argon pour les 

conditions suivantes: p = 1 torr, V = 300 v , f—10 Mhz, et d =2 
cm., a) simulation de M.C, b) ) modèle fluide. 


III. CONCLUSION 


Dans l’ensemble , les prédictions faites par les 
deux modèles (modèle fluide et simulation de Monte 
Carlo) , concernant la densité électronique , le champ 
électrique , le courant de déplacement , le courant 
électronique et le courant ionique sont 
qualitativement et quantitativement en bon accord . 
Ce qui prouve la validité du modèle fluide , utilisé 
dans ce travail , à la pression de 1 torr . 
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